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Abstract. We have searched for the pair production of first generation scalar 
leptoquarks using the full data set (123 pb _1 ) collected with the D0 detector at 
the Fermilab Tevatron during 1992-1996. We observe no candidates, consistent 
with the expected background. We combine these new results from the ee + jets 
and ev + jets channels with the published vv + jets result to obtain a 95% CL 
upper limits on the LQ pair production cross section as a function of mass and 
(3, the branching ratio to a charged lepton and a quark. Comparing to the NLO 
theory predictions, we set 95% CL lower limits on the LQ mass of 225, 204, and 
79 GeV/c 2 for j3 = 1, ^, and 0, respectively. The results of this analysis rule out 
an interpretation of the excess of high Q 2 events at HERA as leptoquarks with LQ 
mass below 200 GeV/c 2 for values of j3 > 0.4. 

1 Introduction 

Leptoquarks (LQ) are hypothesized exotic color-triplet bosons which couple 
to both quarks and leptons. They appear in extended gauge theories and 
composite models and have attributes of both quarks and leptons such as 
color, fractional electric charge, and lepton and baryon quantum numbers 
[0. Leptoquarks with universal couplings to all flavors would give rise to 
flavor-changing neutral currents and are severely constrained by studies of 
low energy phenomena Therefore, only leptoquarks which couple within 
a single generation are considered. The HI and ZEUS experiments at HERA 
have reported an excess of events at high Q 2 in e + p collisions H Q. One 
possible interpretation of these events is production of first generation lepto- 
quarks at a mass near 200 GeV/c 2 || . 

Leptoquarks would be dominantly pair-produced via strong interactions 
in pp collisions, independently of the unknown LQ-/-q Yukawa coupling. 
Each leptoquark would subsequently decay to a lepton and a quark. For 
first generation leptoquarks, this leads to three possible final states: ee+jets, 
ev+jets and vv+jets, with rates proportional to /3 2 , 2/3(1 — (3) and (1 — (3) 2 , 
respectively, where (3 denotes the branching fraction of a leptoquark to an 
electron and a quark (jet). 

Lower limits on the mass of a first generation leptoquark were published 
by the LEP experiments §, by CDF and D0 §, and by HI f§. 

1 Presented at the International Europhysics Conference on High Energy 
Physics (EPS'97), Jerusalem, Israel, Aug. 19-26, 1997. 
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This report describes a search for the pair production of first generation 
scalar leptoquarks in the ee 4- jets and ev + jets final states using 123 ± 7 
pb _1 of data collected by D0 at the Fermilab Tevatron with y/s = 1.8 TeV 
during 1992-1996. A similar search was conducted by the CDF collaboration 
[ 3L The D0 results reported in this article are described in Rcfs. [[it)) and 
[Ly. The D0 detector and data acquisition system are described in detail in 
Ref. @. 

2 The ee + jets Channel 

A base data sample of 101 events with two electrons and two or more jets 
was selected. The electrons (E T > 20 GeV) were required to be separated 
from jets (E T > 15 GeV). Events whose ee invariant mass lies within the Z 
boson mass region were rejected. The efficiency of the trigger used to collect 
the base data sample exceeded 99% for the leptoquark mass range addressed 
by this analysis. 

Monte Carlo (MC) signal samples were generated for leptoquark masses 
between 120 and 260 GeV/c 2 using the ISAJET event generator and a detec- 
tor simulation based on the GEANT program. Leptoquark production cross 
sections were taken from the recently available next-to-leading order (NLO) 
calculations [ [l3| . The primary backgrounds to the ee + jets decay mode are 
Drell-Yan+2 jets production (DY), tt production, and misidentified multijet 
events. The ISAJET DY cross section normalization was fixed by comparing 
it with Z + 2 jets data. The D0 measured tt production cross section of 
5.5 ± 1.8 pb was used for top quark MC HERWIG events. The multijet 
background was estimated directly from data (l(| . 

To search for leptoquarks, a random grid search method was used to 
optimize cuts on the data and MC samples. Consistent results were obtained 
using a neural network |jq ] . The limit setting criterion of a maximum number 
of signal events for a fixed number of background events was adopted. The 
background level chosen was 0.4 events, corresponding to a 67% probability 
that no such events would be observed. 

The set of cuts which optimally separates signal from background was 
determined by a systematic search using distributions of signal MC events. 
Many sets of selection criteria were explored. A cut on a single, relatively 
simple variable, St = H T + H 3 T , where H T = E T X + E T 2 and H T = J2j C ts ^ti 
satisified the limit setting criterion. Approximately 0.4 background events 
are expected for St > 350 GeV. No events remain in the base data sample 
after this St cut is applied. The highest value of St seen in the data is 312 
GeV. The background which is roughly equally distributed among the three 
main sources is estimated to be 0.44 ± 0.06 events. 

To investigate the background further, constrained mass fits were per- 
formed on the events in the base data sample, on background samples, and on 
the 200 GeV/c 2 leptoquark signal MC sample. The 3C mass fit was based on 
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Fig. 1. St vs. 3C fit mass distributions for (a) background, (b) 200 GeV/c 2 lepto- 
quarks, and (c) the base data sample, (d) Mass distribution of the data events (solid 
circles), expected background (solid line histogram), and 200 GeV/c 2 leptoquarks 
(hatched histogram). The inset plot shows these distributions for St > 250 GeV. 



the SQUAW kinematic mass fitting program and required the two ej masses to 
be identical. Figures l(a-c) show St as a function of the fit mass for the esti- 
mated background, 200 GeV/c 2 leptoquark events, and the base data sample. 
The distribution from the data agrees with that of the expected background. 
Figure 1(d) shows the one dimensional mass distributions for the same sam- 
ples. Inset in Fig. 1(d) are the same distributions after a cut on St > 250 
GeV. As can be seen, the data are consistent with the background prediction. 

The overall signal detection efficiency is 9-37% for leptoquark masses 
of 120-250 GeV/c 2 . A 95% confidence level (CL) upper limit on the cross 
section a was set using a Bayesian approach with a flat prior distribution 
for the signal cross section. The statistical and systematic uncertainties on 
the efficiency, the integrated luminosity, and the background estimation were 
included in the limit calculation with Gaussian prior distributions. The re- 
sulting upper limit on the cross section is shown in Fig. 2 together with the 
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NLO calculation of Rcf. Jll|. The intersection of our limit curve with the 
lower edge of the theory band (fx — 2Mlq) is at a = 0.068 pb, leading to a 
lower limit on the mass of a first generation scalar leptoquark of 225 GeV/c 2 . 
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Fig. 2. Upper limit on the leptoquark pair production cross section for 100% decay 
to eq. Also shown is the NLO calculation where the central solid line corre- 
sponds to n = Mlq, and the lower and upper dashed lines to pi = 2Mlq and 
/j, = Mlq/2, respectively. 



3 The eu + jets Channel 

A base data sample of 14 events with one electron, two or more jets, missing 
E T (tf T > 30 GeV), and M T V > 110 GeV/c 2 was selected. The estimated 
background is 17.8 ±2.1 events, of which 11.7 ± 1.8, 4.1 ± 0.9, and 2.0 ±0.7 
events are from W + jets, QCD multijets, and tt production, respectively. 

Monte Carlo samples similar to those in the ee + jets were used in this 
analysis. The W + jets background which is dominant prior to requiring 
M T V > 110 GeV/c 2 was simulated using the VECBOS program. 

Two variables that provide significant discrimination between signal and 
remaining background were identified. They are the energy variable St = 
E T + E T X + E J T 2 + I$t , where E T X and E 3 T 2 are the transverse energies of the 

two jets, and a mass variable ^j-(-Mlq) = min(^^g-^^;i = 1,2), where 

Mlq is an assumed leptoquark mass and M e ( 9 are the invariant masses of 
the electron with the two jets. 
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To find the optimal selection cuts, the same criterion as in the ee + jets 
analysis jio) of maximizing the signal efficiency for a fixed background of «0.4 
events was adopted. In the low mass range (Mlq <120 GeV/c 2 ), where the 
LQ production rates are high, requiring Sj- > 400 GeV is sufficient. For 
Mlq > 120 GeV/c 2 , neural networks (NN) were used since they provide 
higher efficiency than an St cut alone. At each mass where MC events were 
generated, a three layer feed-forward neural network was used with two 
inputs (St and ^^(Mlq)), five hidden nodes, and one output node. Each NN 
was trained using simulated LQ events as the signal (desired output T>n n = 
1) and a mixture of tt , W + jets, and multijet events as background (with 
desired output T>nn = 0). Cuts on T>nn that yield background estimates 
in closest proximity to the desired background were obtained by varying 
T>nn in steps of 0.05. The expected background after the cut ranges between 
0.29 ± 0.25 and 0.61 ± 0.27. No data events pass the cuts. 

Figures 3 (a)-(c) show the 2-dim. distributions of ^-(180) vs. St for sim- 
ulated LQ signal events with Mlq = 180 GeV/c 2 , the combined background, 
and data. The contours corresponding to constant values oIT>mn demonstrate 
the separation achieved between signal and background. The distribution of 
2? aw for data is compared with the predicted distributions for background 
and signal in Fig. 3 (d). The data are described well by background alone. 
The highest T>mn observed in the final data sample is 0.79. 

Using Bayesian statistics, a 95% CL upper limit on the leptoquark pair 
production cross section was obtained for (3 = \ as a function of leptoquark 
mass. The statistical and systematic uncertainties in the efficiency, the in- 
tegrated luminosity, and background estimation were included in the limit 
calculation with Gaussian prior probabilities. The measured 95% CL cross 
section upper limits for various LQ masses are also plotted in Fig. 4 together 
with the NLO calculations (l3) for (3 = | . The intersection of the limit curve 
with the lower edge of the theory band (fi = 2Mlq) is at 0.19 pb, leading to 
a 95% CL lower limit on the LQ mass of 175 GeV/c 2 . 

4 Combined Limits 

An analysis of the vv + jets channel is accomplished by making use of our 
published search (J Ldt m7 A pb^ 1 ) for the supersymmetric partner of the top 
quark (lq . Three events survive the selection criteria (I$t >40 GeV and 2 
jets with Et >30 GeV) consistent with the estimated background of 3.5±1.2 
events. The signal efficiency for Mlq— 80 GeV/c 2 is calculated to be 2.2%. 
This analysis yields the limit Mlq >79 GeV/c 2 at 95% CL for /? =0. 

Combining the ee+jets, e^+jets, and w+jets channels, 95% CL upper 
limits on the LQ pair production cross section were calculated as a function of 
LQ mass for various values of (i. These cross section limits for (3 — h (shown 
in Fig. 4), when compared with NLO theory ([i = 2Mlq), yield a 95% CL 
lower limit on the LQ mass of 204 GeV/c 2 . The lower limits on the LQ mass 
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Fig. 3. ^-(180) vs. St distributions for (a) predicted background, (b) simulated 
LQ events {M L q = 180 GeV/c 2 ), and (c) data, after all cuts except that on 
©jvjv(180). The contours correspond to Vnn = 0.75, 0.85, and 0.95. The area 
of a box is proportional to the number of events in the bin, with the total number 
of events normalized to 115 pb _1 . (d) Distributions of X>jvjv for data (solid circles), 
background (solid hist.) and expected LQ signal for Mlq = 180 GeV/c 2 (hatched 
hist.). 



derived as a function of /?, from all three channels combined as well as from 
the individual channels are shown in Fig. 5. These results can also be used 
to set limits on pair production of any heavy scalar particle decaying into a 
lepton and a quark, in a variety of models. 



5 Conclusion 



We have presented a search for first generation scalar leptoquark pairs in the 
ee+jets and e^+jets decay channels. Combining the results with those from 
the vv + jets channel, we exclude at 95% CL leptoquarks with mass below 
225 GeV/c 2 for = 1, 204 GeV/c 2 for = ±, and 79 GeV/c 2 for = 0. 
Our results exclude (at 95% CL) an interpretation of the HERA high Q 2 
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Fig. 4. Measured 95% CL upper limits on the leptoquark pair production cross 
section in the eu+jets channel (circles) and all three channels combined (triangles) 
for /? = i. Also shown are the NLO calculations of Ref. E3| where the central 
line corresponds to jj, = Mlq, and the lower and upper lines to n — 2Mlq and 
H — Mlq/2, respectively. 

event excess via s-channel LQ production with LQ mass below 200 GeV/c 2 
for values of (3 > 0.4. 
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are also shown. The shaded area is excluded at 95% CL. 
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